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Introduction 

In  Z-pinch  experiments  a  strongly  radi¬ 
ating  plasma  may  adversely  effect  the  trans¬ 
port  of  power  delivered  to  it.  Plasma  gene¬ 
rated  x-rays  cause  ablation  of  the  diode 
walls  and,  if  the  walls  do  not  become  good 
conductors  or  if  the  current  (I)  carried  in 
the  electrode  surfaces  is  too  small,  loss  of 
insulation  in  the  diode  can  result.  Last 
year  we  interpreted  experimental  results  at 
Maxwell  Laboratories!  as  evidence  for  hydro- 
dynamic  closure  and  proposed  a  theory  to  ex¬ 
plain  the  reported  behavior.  While  the 
explanation  was  qualitatively  correct  there 
were  no  measurements  that  could  be  compared 
to  a  detailed  quantitative  analysis.  The 
work  described  here  presents  the  theoretical 
model  used  to  design,  and  be  compared  with, 
recent  experiments  on  BLACKJACK  III  (BJ3)  at 
Maxwell  Laboratories. 

These  experiments  are  required  to  vali¬ 
date  the  theory  so  that  the  effects  of  closure 
can  be  predicted  for  present  and  future  SHIVA 
configurations,  diodes,  and  plasma  loads. 

Analysis  of  BLACKJACK  III  Experiments 
Experimental  Geometry 

Figure  1  shows  a  cylindrical  plasma  ir¬ 
radiating  the  surfaces  which  carry  current  to 
and  from  it.  This  double-disk  geometry  is  the 
simplest  one  in  which  to  analyze  the  radia- 
tively  driven  ablation  and  during  the  design 
phase  for  these  experiments  it  was  selected 
as  providing  a  reasonable  high  power  inductive 
diode  while  simultaneously  providing  a  con¬ 
figuration  in  which  (a)  the  radiative  trans¬ 
port  was  not  too  complicated  and  (b)  the  blow- 
off  from  the  electrodes  could  be  seen  and 
measured  easily. 
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Figure  1.  Double-disk  geometry  schematic. 

The  experimental  configuration  is  shown 
in  Figure  2.  The  load  for  the  diode  is  an 
array  of  thin  wires  which,  on  imploding,  emit 
an  intense  burst  of  x-rays.  The  diode  itself 
consists  of  two  parallel,  flat  disks  (alumi¬ 
num  or  copper)  approximately  25  cm  in  diameter 


With  this  configuration  it  is  possible  to 
have  unobstructed  view  of  the  plasma  as  well 
as  the  movement  of  the  electrode  surfaces. 


Visual  diagnostics  include  streak  photo¬ 
graphy  as  well  as  laser  shadowgraphy .  The 
streak  photograph  information  turns  out  to 
be  very  difficult  to  interpret  and,  in  general 
not  useful  for  this  experiment.  The  laser 
shadowgraph  information,  on  the  other  hand, 
gives  information  about  the  lowest  density 
parts  of  the  ablating  plasma  which  is  useful 
for  comparison  purposes. 

Computational  Model 

The  model  used  to  describe  the  inter¬ 
action  of  the  x-rays  with  the  walls  has  to 
include  the  following  features:  (a)  x-ray 
transport  in  the  vacuum  region;  (b)  energy 
diffusion  into  the  walls;  (c)  opacity  and 
equation  of  state  information  for  the  wall 
material;  (d)  hydrodynamics  of  the  walls; 

(e)  magnetic  forces  on  the  expanding  walls. 

To  do  this  we  used  a  version  of  the  MRVERA 
code2  which  was  modified  to  include  the  mag¬ 
netic  field  terms  (JxB  force,  joule  heating, 
Ohm's  Law,  Faraday's  Law,  Ampere's  Law)  and 
thermal  conduction.  The  individual  sections 
of  the  wall  were  treated  locally  as  one¬ 
dimensional  but  they  were  coupled  through  the 
two-dimensional  features  of  the  vacuum  radia¬ 
tion  transport.  Extensive  use  was  made  of 
the  CAP03  generated  equation  of  state  and 
opacity  data  for  the  electrode  materials 
(aluminum,  plastic  and  copper). 

Preliminary  Calculations 

The  x-ray  output  and  diode  wall  current 
histories  are  shown  in  Figure  3.  In  the  ini¬ 
tial  calculations  the  area  under  the  power 
curve  was  8  kJ  (this  turned  out  to  be  too 

large)  and  I  was  0.6  MA.  The  gap  separa- 
max 

ting  the  disks  was  2  cm. 

Figure  4  shows  the  diode  geometry  as 
well  as  the  calculational  regions'  surface 
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Figure  3.  Normalized  XRD  and  current  pulse- 
iorms . 
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Figure  4.  Plasma  surface  temperature  (in  eV) 
at  115,  200,  and  370  ns,  respec¬ 
tively  . 

temperature  (0)  at  three  times.  The  rela¬ 
tively  cool  temperatures  (“v5  eV)  result  in 
the  emission  of  very  soft  x-rays  in  the  10 
to  20  eV  band  with  an  appreciable  fraction 
radiated  in  the  visible. 

An  interesting  feature  which  resulted 
from  these  calculations  is  the  ablation 
front's  non-monotonic  dropoff  with  radius  (r), 
shown  in  Figure  5.  Early  in  the  pulse,  be¬ 
fore  the  magnetic  field  (B)  builds  up  and  its 
2 

1/r  dependence  dominates,  blow-off  is  most 
rapid  close  to  the  source.  Later  the  B  field 
increase  is  enhanced  by  the  hotter  plasma 
there  being  a  good  conductor,  thus  preventing 
magnetic  diffusion  from  taking  place  and  pro¬ 
ducing  steep  gradients  in  the  magnetic  force, 
o 

3/9z  (B  /8ir)  .  At  late  times  B  has  dropped  off 
and  ceases  to  restrain  the  plasma.  Farther 
from  the  source  the  weaker  B,  coupled  with  a 
cooler  more  resistive  plasma,  has  less  and 
less  effect  on  the  hydrodynamic  motion.  For 

2 

this  geometry  the  flux,  (energy/cm  /sec) 

3 

normal  to  the  surface  drops  off  as  v(i/r) 
for  large  r.  Thus  at  larger  r  the  material  0 
and  blowoff  velocity  (vf)  are  reduced. 

A  simple  quasi-analytic  model  described 
in  a  recent  report^  can  be  used  to  further 


Figure  5.  Blow-off  front  along  2  cm  gap  BJ3 
disk  diode. 

understand  the  phenomenology.  If  the  mag¬ 
netic  field  pressure  can  be  neglected  (i.e., 
far  from  the  source),  the  velocity  (v^)  of 

the  ablation  front  exposed  to  an  incident 
radiant  flux,  dE/dt,  can  be  parameterized  as 

2 

vf  =  vf(E),  where  E  is  the  energy  (J/cm  ) 

that  has  entered  the  material  after  front 
surface  vaporization  has  occurred.  Separate 
radiation  hydrodynamic  calculations  yielded 

2/9 

the  dependence  v^  2.58  E  '  cm/ysec,  which 
holds  (in  AH)  over  a  wide  range  of  conditions. 

To  compute  vf  in  a  simple  manner  the 

regions  are  assumed  to  receive  energy  directly 
from,  and  only  from,  the  source.  Neglect  of 
re-emission  may  be  quite  a  good  approximation 
for  regions  that  are  cold  compared  to  the 
source.  Figure  5  shows  the  analytic  blow-off 
front  profile  histories  compared  to  those 
calculated  with  a  magnetic  field  and  re¬ 
emission  by  MRVERA.  Beyond  r  =  3  cm  the 
analytic  treatment  performs  well  showing  good 
agreement  during  the  MRVERA  calculational 
time  span  for  the  far  field  regions  where  the 
influence  of  B  was  expected  to  be  slight. 

Interpretation  of  Visible  Streak  Data 

The  visible  streak  photographs  give  clear 
evidence  for  the  expansion  of  both  cathode  and 
anode  surfaces  and  for  the  closing  of  the  di¬ 
ode  gap.  Figure  6  shows  the  displacement  of 
the  ablation  front  with  respect  to  the  orig¬ 
inal  cathode  and  anode  surfaces  observed  from 
a  visible  streak  on  B J3 .  The  difference  in 
the  cathode  and  anode  blow-off  rates  are 
slight,  probably  due  to  misalignment  and  not 
considered  to  be  experimentally  significant. 

It  is  clear  that  the  expansion  rate  is  fairly 
constant  between  100  and  200  nanoseconds  (ns) . 
Thereafter  it  slows  down  or  comes  to  a  halt . 
This  behavior  is  not  expected  theoretically 
(Figure  5). 

The  reason  for  this  behavior,  it  turns 
out,  is  based  on  the  interpretat ion  of  what 
is  observed.  Briefly,  the  visible  streak 
photographs  do  not  provide  direct  information 
about  the  displacement  of  the  electrode  sur¬ 
faces.  Rather,  as  we  have  shown  elsewhere, 
they  provide  information  about  the  more 
strongly  emitting  regions  of  the  plasma. 

Beyond  this,  considerable  unfolding  is  re¬ 
quired  to  obtain  anything  that  is  useful  for 
comparison  purposes. 
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Figure  6.  Displacements  observed  in  BJ3 
along  the  chord  y  =  1.5  cm 
(visible  streak  photographs). 


Figure  7.  Comparison  of  theory  and  experi¬ 
ment  for  the  maximum  ablation 
front  displacement  of  electrode 
surfaces  in  BJ3  (laser  shadow¬ 
graph)  . 


analyze  these  since  it  is  a  major  task  to 
unfold  the  measurements  and  relate  them  to 
calculated  quantities.  However,  those  pro¬ 
files  do  exhibit  the  non-monotonic  dropoff 
of  displacement  amplitude  z(r)  with  radius. 


Interpretation  of  Laser  Shadowgraph  Results 

The  other  experimental  diagnostic  which 
was  used  in  the  BLACKJACK  III  experiments  was 
a  laser  shadowgraph.  In  this  technique  a 
visible  laser  provides  the  source  to  be  meas¬ 
ured;  whereas,  the  visible  streak  relied  on 
material  emission  sources.  The  advantage  of 
this  scheme  is  that  the  top  of  the  ablated 
layer  can  now  be  seen.  Even  at  densities 
where  the  material  is  no  longer  opaque  it  can 
be  "seen"  because  the  medium  refracts  the 
visible  light  passing  through  it.  The  light 
that  is  strongly  refracted  due  to  large  elec¬ 
tron  density  (ng)  gradients  will  not  show  up 

on  the  receiving  system.  The  technique  high¬ 
lights  strong  ng  gradients  and  is  thus  re¬ 
liable  for  observing  the  ablation  front  if  it 
is  not  too  diffuse. 

The  shadowgraph  method  was  used  with  the 
double-disk  geometry  with  radial  viewing,  as 
shown  in  Figure  2,  and  shadowgraphs  were  snap¬ 
shot  at  intervals  of  30-50  ns.  It  was  used 
with  ASL ,  brass,  and  Cu  electrodes  and  almost 
entirely  with  diode  gaps  of  1  cm  at  atmo¬ 
spheric  pressure  and  8  mm  under  vacuum  con¬ 
ditions  . 

MRVERA  performed  calculations  using  Cu 
electrodes  and  an  8  mm  gap.  The  energy  under 
the  power  curve  (Figure  3)  was  the  measured 
value  of  2.5  kj.  Figure  7  shows  a  comparison 
of  measurements  on  Cu,  brass,  and  AS,  with  the 
calculated  results  for  Cu.  The  maximum  dis¬ 
placement  (Az)  of  the  anode  and  cathode  sur¬ 
faces  are  plotted  as  a  function  of  time. 

These  points  are  special  in  the  sense  that 
they  are  unobscured  by  material  in  the  line 
of  sight.  The  radial  positions  corresponding 
to  the  points  indicated  are,  in  general,  time- 
dependent.  This  is  to  be  expected  since,  in 
Figure  5,  Az  occurs  at  larger  r  late  in  the 

Hid,  A 

pulse.  Agreement  between  theory  and  experi¬ 
ment  is  remarkably  good,  considering  the  com¬ 
plexities  of  the  underlying  physics. 

Not  shown  here  are  the  projected  radial 
profiles  of  the  shadowgraph.  We  did  not 


It  is  interesting  that  the  cathode  and 
anode  blow-off  are  not  quite  symmetric.  Our 
theory  says  that  if  the  plasma  source  is 
uniform  across  the  gap  then  both  should  show 
equal  expansion.  A  possible  explanation  for 
the  slight  asymmetry  is  the  difficulty  which 
is  intrinsic  to  aligning  the  laser  apparatus 
and  the  beam  parallel  to  the  electrode  sur¬ 
face.  Deviations  from  perfect  alignment 
could  result  in  the  observed  results.  The 
angle  subtended  by  the  blow-off  is  40  ^  2  mm/ 

_3 

100  cm  ^  2  x  io  radians  or  a  tenth  of  a 
degree.  The  asymmetry  errors  involve  even 
smaller  angles. 

The  other  feature  which  shows  deviation 
from  theory  occurs  late  in  the  pulse.  Theory 
predicts  closure;  whereas,  the  experimental 
results  show  a  slowing  of  the  closure  rate. 

We  suspect  the  reason  for  this  is  again  the 
difficulties  of  the  experiment.  Until  150 

or  200  ns  the  B  field  ensures  that  the  n  is 

e 

very  steep  near  the  blow-off  front  (Figure  8) 
so  that  the  shadowgraph  technique  is  ideal. 


DISTANCE  (MM) 

Figure  8.  Electron  density  profiles  at  the 
radial  positions  indicated  both 
early  and  late  in  the  pulse. 
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Later,  I  and  B  drop  and  the  front  expands  and 
becomes  quite  diffuse.  The  reduced  ng  gradi¬ 
ents  are  much  more  difficult  to  measure  and, 
we  suspect,  the  refraction  too  weak  to  indi¬ 
cate  the  low  n  front  location, 
e 

Agreement  between  calculations  and 
measurements  for  Cu  is  excellent.  However, 
measurements  for  A&  and  brass  electrodes  are 
not  markedly  different.  Reasons  for  this 
behavior  are  discussed  elsewhere. Ref¬ 
erence  6  also  examines  the  effects  of  surface 
contaminants . 

Conclusions 

We  have  demonstrated  that  for  high  power 
diodes  in  which  the  load  is  not  passive  but 
emits  x-rays,  the  wall  electrodes  undergo  ex¬ 
tensive  ablation.  The  effect  is  thus  impor¬ 
tant  in  all  high  power  Z-pinch  experiments  in 
plasmas.  We  have  seen  how  closure  could  take 
place  and,  while  in  these  experiments  complete 
closure  did  not  occur  until  late  in  the  pulse 
due  to  a  careful  choosing  of  the  pulse  length 
and  gap  width,  it  is  clear  that  without  spe¬ 
cial  precautions  closure  can  occur  during  ex¬ 
periments. 

Even  without  the  complete  closure  of  the 
diode  the  power  transfer  to  the  load  can  be 
affected  appreciably.  That  is,  the  local  in¬ 
ductance  (per  unit  length)  or,  equivalently, 
the  local  line  impedance,  is  altered  by  the 
moving  electrode  surfaces.  This,  in  turn, 
modifies  the  voltage  and  hence  the  current 
carried.  The  effect  is  discussed  elsewhere.8 

The  computational  model  appears  to  have 
been  verified  by  experiment.  This  is  an  im¬ 
pressive  achievement  in  a  physical  regime 
where  the  transport  coefficients,  opacities, 
etc.  are  somewhat  suspect.  Moreover,  the  re¬ 
sults  obtained  show  that  the  radiation-MHD 
techniques  used  to  compute  fluid  behavior  in 
this  very  high  aspect  ratio  geometry  are  valid. 
This  provides  encouragement  for  theoretically 
predicting  the  closure  properties  of  x-ray 
diodes  and  for  designing  them  for  specified 
radiating  loads.  It  also  provides  incentive 
for  checking  the  thermodynamic,  electrical, 
and  opacity  behavior  of  materials  other  than 
those  considered  here.  Finally,  the  results 
show  that  the  laser  shadowgraph  (or  holograph) 
technique  is  invaluable  for  observing  these 
cool  diffuse  plasmas  and  is  superior  to  visi¬ 
ble  streak  diagnostic  photographs. 
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